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We give a method for determining the distortion effect of a nonlinearity 
in a feedback loop. 

I. INTRODUCTION 

Desoer gives an interesting analysis of distortion resulting from 
a nonlinearity of the form v + ev m {m an odd integer) in a feedback 
loop. 1 Sandberg considers virtually the same problem for nonlineari- 
ties with upper and lower bounds on the slope. 2 On page 2546 of his 
work, Sandberg suggests that Desoer's result may be sharpened. 
Our purpose is to show how a small modification of Desoer's analysis 
might give this sharpening and extend its applicability. 

Desoer's method is to find conditions for a particular mapping to 
be a contraction in a ball. The method presented in another work 
is particularly suited to that problem and will be applied in this 
paper. 3 The problem of distortion in nonlinear systems is also con- 
sidered in References 4 and 5 among other papers. 

II. NOTATION AND PRELIMINARIES 

The feedback loop (with unity feedback for simplicity) is assumed 
to be described by 

y = NL(r - y) (1) 

where the input r and output y arc in some Banach space. L and N 
are linear and nonlinear operators, respectively, mapping the Banach 
space into itself. We need not, at this point, specify which Banach 
space we are working in. Rather, we refer the reader to Reference 2 
for details on two Banach spaces of interest for analysis of nonlinear 
feedback loops." In particular, Reference 2 shows how to evaluate 

*It must bo verified that the Banach space (or an appropriate subset) is 
mapped into itself by the nonlinearity. In particular, nonlinearities such as de- 
scribed by polynomials do not map La into itself. 
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the norm of the linear operator when it is denned by a convolution 
operation or by a transfer function (frequency response). 

in. THE PROBLEM OF DISTORTION 

Suppose that 

N(x) =x+ eP(x). (2) 

Then the loop is linear if e = and it is of interest to determine how 
the loop response differs for a nonzero e. This difference is called the 
distortion. On the other hand, we might consider some fixed | c | > 
and determine how small the input r must be in order that the dis- 
tortion is sufficiently small. This latter question assumes that P{x) 
is of an order less than x as x— »0. 

The following manipulation is convenient for this problem. From 
equations 1 and 2 we have 

y = L(r - y) + eP[L(r - y)]. (3) 

If we assume that (I+L) has a bounded inverse where I is the identity 
map,* we obtain 

y = (I + L)~ l Lr + (I + L)" l eP[L(r - y)\. 
Then, if z is the response of the linearized loop, 

z = (I + LY'Lr. 
And if £ represents the distortion, 

$ - y - z, 

we have 

{ = 6(1 + Ly'Plz - m 

We are thus interested in finding a fixed point of the operation M (£) . 
In particular, how large is £? To solve this problem, we use a con- 
venient modification of the contraction mapping fixed point theorem. 

IV. THE CONTRACTION MAPPING THEOREM 

Let X be a complete metric space (with metric d) containing 
the closed set Q and let F map £2 into itself. F is a contraction map- 
ping if there is an a t [0, 1) such that 

d[F(x), F(x')] ^ ad(x, x') (x, x' e Q). 

♦For conditions for the existence of this bounded inverse, see Reference 2, 
especially p. 2538, 
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The contraction mapping theorem (Reference 6, p. 627) states that 
if F is a contraction mapping then there is a unique x* £ fi such that 
x* = F(x*), that is, x* is a fixed point of the operation F. Also, x* 
is the limit of a sequence {x n } where 

.r n + l = F(x n ) 
and x is any element of a. 

One aspect of using the above theorem is finding the appropriate 
set n mapped into itself. Often, the contraction mapping theorem is 
used when n is the whole space, that is, F is globally Lipschitzian. 
The analysis of Reference 1 may be viewed as a method of deter- 
mined a ball about the origin such that a mapping is a contraction 
in that ball. The general problem of simultaneously trying to de- 
termine a set mapped into itself such that a mapping is contraction 
on that set is discussed in Reference 3. The following simple theorem 
from Reference 3 is useful in this direction. 

Theorem: Let B be a Banach space. F maps B into itself and x Q t B. 
It is assumed that 

(i) F has a derivative at all x z B 
(it) There is a nondecreasing function g such that if x e B, then 

II F'(x) |! ^ 0(|| x - Xo ||) 

(Hi) There is an a e [0, 1) such that 

k 



1 — a/ 

where 

k^ \\F(x ) -x \\. 
Then there is a unique x* 1 12 such that 

x* = F(x*) 
where 

n = l [ x:xzB, || x- x.W^^-A. 

Remarks: See chapter XVII of Reference 6 for a general discussion 
of differentiation in Banach spaces. 

It is often a straightforward matter to find an appropriate func- 
tion g as we shall see in the distortion problem of this paper. 
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V. SOLUTION OF THE DISTORTION PROBLEM 

To apply the preceding theorem to the distortion problem of Sec- 
tion III, we first find M'{$), then a nondecreasing g such that 

\\M'®\\ Sf(iie-«.||) -rfll ill) fo-o). (4) 

And with P (0) =0 (for simplicity) , 

II M%) - 5o I! = II e(/ + L)"P( 2 ) || 

g fc. (5) 

We must finally find an a e[0.1) such that 

'' ^ «. (6) 



\\ - a 

With (7+L)- 1 and L both assumed to be bounded linear operators, 
we have 

II Jf'ft) II = II e(7 + Ly l P'(z - WL l| (7) 

(assuming that P has a Frechet derivative) . It should be clear that 
our method of analysis is not restricted to nonlinearities described 
by functions of the form of v + ev m as used in Reference 1. For the 
case of the space of continuous real valued functions with the sup 
norm* and 

p(x) = x m m an integer > 1 (not necessarily odd) (8) 

we have that 

P'(x) = mx m ~\ (9) 

Then, using equations 7 and 9, 

II Jf'ft) II £ I « I -II (I + L)- 1 ||-m-|| (z - W'- 1 \\-\\L || 

^ W |e hi (/ + £>-' ii (ii « n + ii i ii-ii siir"* 1 II ^11 

- 9(11 € ~ So II) 

= 0(11 € II) Go - 0). (10) 

Now, using equations 5 and 6, we obtain 

» | • |-|| (J +.13"' II 

, l|L|l (i MI + i i Mi-NN i a_^)-i i-i ^ i r )"- ga , (11) 

* What might be considered to be a disadvantage of using this space is that 
the norms of the linear operator are expressed in terms of impulse responses 
rather than frequency responses. 
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Condition (11) could be used in several ways. For a fixed e, we could 
determine how small || z || has to be in order for there to be an 
a t [0, 1) satisfying (11) and thus get a bound on the distortion || £ ||. 
The emphasis in Reference 1 is in determining how small e should be 
with linearized outputs satisfying \\z\\ < 1 and the distortion \\ £ \\ < % 
in order for the method of successive approximations to converge at 
a given rate (a = j). The discussion on page 2546 of Sandberg's article 2 
assumes the following conditions (in our notation) : 

m = 3 

II L || = 100 (12) 

|| (/ 4- L)- 1 || = 2. 

Then, (11) becomes 

600 | . |(l + m^J S i. (13) 

If |e| is less than about 1/2900 (actually a little larger, then (13) 
is satisfied. Then, the distortion | satisfies 

II til* l 



1 -a 

^ leUllcr + LrMMMp 

stM. (14) 

The bound obtained using equation 25 of Desoer's article 1 is | e | ^ 
1/(2150-2900), a substantially smaller bound. 

VI. CONCLUSION 

Notice that since we do not require the mapping to be a contrac- 
tion in the whole space, we only get uniqueness in Q, the ball of 
radius k/(l—a). However, the result may be strengthened by also 
seeking the largest contraction constant a, satisfying condition in of 
the theorem. Then the fixed point is also unique in the larger ball. 
On the other hand, uniqueness information might be available from 
another source (for example, a property of a differential equation). 

We notice that the existence of derivatives in the theorem may 
actually be relaxed if there is a nondecreasing function suitably 
bounding Lipschitz constants. We also mention the possibility of using 
transformations to facilitate the application of the result. 
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The following may be helpful in visualizing the application of the 
contraction mapping theorem.* Assume that condition Hi of the 
theorem of Section IV is satisfied with equality, that is, 



k 



= a. 



a - a, 

The radius of the ball Q, is fc/(l - a). Letting 

k 



r = 



1 - 



the condition is seen to be 



g(r) = i - - 



which Fig. 1 shows pictorially. 




-||P(Xo)-«oH 



RADIUS OF BALL CONTAINING 

ALL ITERATES Xl +1 = F(I L ) 

AND A FIXED POINT 



g(r): LIPSCHITZ 
CONSTANT ON BALL 
OF RADIUS r 



r= x-io 



c— BALL OF GUARANTEED UNIQUENESS— >j 

Fig. 1 — Contraction mapping theorem. 
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